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ABSTRACT 



Electric discharges were detected directly in the cloudy atmospheres of Earth, 
Jupiter and Saturn, are debatable for Venus, and indirectly inferred for Neptune 
and Uranus in our solar system. Sprites (and other types of transient luminous 
events) have been detected only on Earth, and are theoretically predicted for 
Jupiter, Saturn and Venus. Cloud formation is a common phenomenon in ultra- 
cool atmospheres such as in Brown Dwarf and extrasolar planetary atmospheres. 
Cloud particles can be expected to carry considerable charges which may trig- 
ger discharge events via small-scale processes between individual cloud particles 
(intra-cloud discharges) or large-scale processes between clouds (inter-cloud dis- 
charges). We investigate electrostatic breakdown characteristics, like critical field 
strengths and critical charge densities per surface, to demonstrate under which 
conditions mineral clouds undergo electric discharge events which may trigger 
or be responsible for sporadic X-ray emission. We apply results from our ki- 
netic dust cloud formation model that is part of the Drift-Phoenix model 
atmosphere simulations. We present a first investigation of the dependence of 
the breakdown conditions in Brown Dwarf and giant gas exoplanets on the local 
gas-phase chemistry, the effective temperature and primordial gas-phase metal- 
licity. Our results suggest that different intra-cloud discharge processes dominate 
at different heights inside mineral clouds: local coronal (point discharges) and 
small-scale sparks at the bottom region of the cloud where the gas density is high, 
and flow discharges and large-scale sparks near, and maybe above, the cloud top. 
The comparison of the thermal degree of ionisation and the number density of 
cloud particles allows us to suggest the efficiency with which discharges will occur 
in planetary atmospheres. 
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1. Introduction 

Electric discharges, including transient luminous event^j] or transient discharges in 
gases, are non-equilibrium plasmas produced by short-lived electrical conditions of which 
lightning and associated sprites are the most widely known. Lightning and sprites are 
commonly observed on Earth. Lightning is also detected or indirectly inferred on most of 
the cloud-carrying solar system planets. Observations on Earth show that cloud particles 
or droplets trigger such transient discharges in atmospheric environments (e.g. Yuan et al. 
2011). It is not obvious that the underlying physical mechanisms for triggering lightning 
and sprites are different: they each involve an unbalanced evolution of accelerated free 
electron populations, driven by electric potential differences. How the latter are created 
from the organisation of ambient gas/particles may be different, but the underlying science 
is not. 

Most of the solar system planets form clouds made of liquids and ices, but dependent 
on local conditions, mineral hazes may also form (e.g. Saunders & Plane 2011). Recent 
observations of HD 189733b (Pont et al. 2008; Sing et al. 2009, 2011; Gibson et al. 2012) 
suggest that extraterrestrial giant gas planets contain small mineral dust particles high up 
in their atmospheres. Further examples include GJ1214b and Kepler 7b. 

Prompted by the abundance of discharge processes in clouds on Earth and in the 
clouds of solar system planets, we take a first look at the conditions for such discharge 
events in extrasolar, planetary and Brown Dwarf atmospheres. The essential requirements 

x In the community of terrestrial atmospheric physics, the term transient luminous event 
is used to describe sprites, elves, and blue-jets only, but not lightning. 
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are the presence of both, an electric field (electric potential difference) and free electrons. 
This electric field accelerates the free electrons through the ambient gas. If they collide with 
other gas particles sufficiently frequently and with enough energy to liberate increasingly 
more electrons, an electron avalanche is produced. The critical electric field required for this 
depends, to some extent, on the composition of the gas. This avalanche can develop into a 
self-propagating ionisation front (streamer) if the self-field of the electrons in the avalanche 
is similar to the strength of the applied external field. If the avalanche has reached a 
critical size, it can propagate into areas with a lower electric field strength in which an 
avalanche would not start normally. Two electrical breakdown mechanisms are known to 
occur in dielectrics: the conventional breakdown occurring for example in glow discharges; 
and the runaway breakdown which is suggested to play a role in lightning discharges on 
Earth. The threshold electric field needed to initiate the avalanche of a runaway breakdown 
appears one order of magnitude below that for a conventional breakdown (for more details 
see Roussel-Dupre et al. 2008; Marshall et al. 1995). Each of the resulting streamers is a 
highly conducting plasma channel in which the density of non-thermal electrons is greatly 
increased, and from which electrons leak into the ambient gas. 

Helling, Jardine & Mokler (2011) have argued that dust clouds that are made of 
mineral particles can be charged e.g. by gas turbulence induced dust-dust collisions. 
The cloud particles remain charged for long enough that discharge events can occur and 
streamers can be established. Other possible processes for charge separation, which we have 
not considered, are e.g. different contact potentials of the material making up the grain, 
and fractoemission or collisional charging of polarised grains in an external electric field 
(Pahtz, Herrman & Shinbrot 2010). The Coulomb recombination time of cloud particles by 
thermal electrons is long enough that electric fields can establish between passing grains 
which then can initiate subsequent electron avalanche processes leading to streamers (Fig. 
4 in Helling et al. 2011). The result is that during the electron avalanche phase there is 
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an exponential increase in free charges for each free electron that exists in the gas phase. 
Numerical experiments that simulate an avalanche from the initial electron up to and 
beyond the avalanche-to-streamer transition demonstrate consistently that a single electron 
seeding a streamer can cause the occurrence of 10 12 — 10 13 secondary non-thermal electrons 
for a certain time interval (Dowds et al. 2003, Ebert et al. 2010). 

Laboratory studies suggest that during grain-grain impact, smaller grains are more 
negatively charged compared to larger grains of the same composition (Lacks & Levandovsky 
2007, Forward et al. 2009, Merrison et al. 2012). A similar conclusion was reached from 
volcano ash experiments (Hatakeyama & Uchikawa 1952, Kikuchi & Endoh 1982, James 
et al. 2000). Experimental results for different minerals are not available at present. 
Experiments on wind-driven entrainment of Martian dust (mainly iron oxides) and Basalt 
dust (Na/K/Ca-silicates) by Merrison et al. (2012) suggest average electric charge values of 
10 3 ... 10 5 e~/grain. The exact values might, however, be influenced by the experimental 
set-up (Aplin et al. 2012). These experiments were carried out under pressures between 
10 mbar and 1 bar. Only a very weak pressure dependence on the number of charges per 
grain was seen. No strong dependence on the grain material was found either, but one 
order of magnitude variation is suggested due to grain sizes, morphology and abundance 
of grains. They found no evidence for a dependence of dust electrification on atmospheric 
composition. 

The extrasolar atmospheric environments we are looking at in Brown Dwarfs and giant 
gas planets have local temperatures between T gas = 500K - 4000K and gas pressures between 
Pgas ~ 10~ 10 bar - 10 bar. This pressure range comprises the pressure interval for Mars-dust 
analog experiments by Merrison et al. (2012). The resulting degree of ionisation by thermal 
processes in Brown Dwarfs and giant gas planet atmospheres is Xc — Pc/Pga,s ~ 10~ 15 - 
10~ 7 throughout the atmosphere, which is very low (Helling et al. 2012), but enough to 
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start an electron avalanche. The formation of dust by seed formation and bulk growth 
takes place in a temperature window of ~500K-~2100K and leads to the formation of 
mineral clouds. Gravitational settling, convective mixing and element depletion are major 
associated processes. Helling, Woitke & Thi (2008) have shown that the upper cloud, (i.e. 
low temperature and low-pressure) will be dominated by small, dirty (i.e. inclusions of other 
materials) silicate grains with inclusions of iron and metal oxides, and the warmer, denser 
cloud base by bigger, dirty iron grains with metal inclusions. The actual size of the cloud 
particles deviate from this mean value according to a height dependent size distribution 
(Fig. 8 in Helling, Woitke & Thi 2008). The chemical composition of the grains as well as 
their size distribution change with height inside a cloud in a quasi-stationary environment. 

The different scale regimes involved in the formation of atmospheric clouds suggest the 
consideration of two different scale regimes with respect to charge and discharge processes. 
The large scale regime is where the cloud particles settle gravitationally (rain out) according 
to their mass and size, or where they are transported by hydrodynamical processes like 
convection. This will leave the smallest and potentially more negatively charged particles 
suspended in the higher layers, while the larger and more positively charged particles 
populate the lower part of the cloud. The small scale regime describes the scales at which 
cloud particles interact during grain-grain collisions, where colliding particles (droplet) or 
individual cloud particles produce point discharges (also called corona) that may start 
large-scale lightning. This analogy is drawn from the discharge experiments for Earth 
conditions where colliding drops in storms and ice particles can produce point discharges 
which are thought to be able to initiate lightning (for more details see MacGorman & Rust 
1998). A similar scenario is suggested in Merrison et al. (2012), and was utilised by Farrell 
et al. (1999) and Michael et al. (2008) to model Martian dust storms. 

What are the conditions to trigger a discharge event which subsequently increases the 
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local density of free electrons by non-thermal processes leading potentially to the occurrence 
of lightning or associated transient luminous events like sprites? The local electric field, E, 
must exceed the ionisation threshold of the ambient gas, hence a critical field strength, E CI it, 
and seed electrons for the development of an electron avalanche need to be present. We note 
that measurements of the actual threshold electric field strength for a terrestrial lightning 
flash in a cloud are not easy to obtain (Marshall et al. 2005). Marshall et al. (1995) 
performed balloon experiments from which they derived electric breakdown thresholds that 
led them to hypothesised that lightning flashes might be initiated by runaway breakdown. 
This paper therefore presents a first investigation of the discharge behaviour in extrasolar 
atmospheres by studying the following questions. 

- What electric field strength is needed to allow breakdown phenomena known from Earth 
and solar system planets to occur in extrasolar, planetary objects and Brown Dwarfs? 

- What can we learn about the breakdown conditions in extrasolar objects, for example 
breakdown distances, minimum voltages, or chemistry dependence? 

- What are the critical conditions for breakdown such as cloud size, grain size, number of 
charges per surface? Under which conditions can the critical field strength be produced by 
grains alone? 

- How many thermal electrons are available in a mineral cloud of a Brown Dwarf or a giant 
gas planet in order to trigger a streamer event? How efficiently can discharge events occur 
and, hence influence the gas phase composition? 

To answer these questions, we first outline the basic ideas and relations for an 
electrostatic breakdown in gases of different chemical compositions in Sect. [2} These 
relations are evaluated in Sect. [3] for previously calculated model atmosphere structures of 
Brown Dwarfs and giant gas planets. We are using Drift-Phoenix model atmosphere 
structures which are the result of the solution of the coupled equations of radiative transfer, 
convective energy transport (modelled by mixing length theory), chemical equilibrium 
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(modelled by laws of mass action), hydrostatic equilibrium, and dust cloud formation (Dehn 
2007, Helling et al. 2008a,b; Witte, Helling & Hauschildt 2009; Witte et al. 2011). The dust 
cloud formation model includes a model for seed formation (nucleation), surface growth and 
evaporation of mixed materials. The effects of gravitational settling (drift) and convective 
overshooting on the cloud formation are explicitly included in our dust model equations 
(Woitke & Helling 2003, 2004; Helling & Woitke 2006; and Helling, Woitke & Thi 2008). 
The results of the Drift-Phoenix model atmosphere simulations include for example 
the gas temperature - pressure structure (T gas , p gas ), the local gas-phase composition, 
the local electron number density (n e ), the number of dust grains (rid) and their height 
dependent sizes (a). These models are determined by the effective temperature, T e g [K], 
the surface gravity, g [cm/s 2 ], and the initial element abundances which are set to the solar 
values unless specified otherwise. The final element abundances are determined by element 
depletion due to dust formation. We also discuss how many charges per surface would be 



needed to initiate a field breakdown based on this model output (Sect. 3.2), and we estimate 



the enrichment rate of the gas-phase by streamer electrons (Sect. 3.4). Section |4j contains 
our conclusions. 



2. Breakdown characteristics 

Electric breakdown is the process of transforming a non-conducting substance into a 
conductor as the result of applying a sufficiently large electric field. This ionisation state 
of the bulk gas phase builds up over times between 10~ 8 . . . 10~ 4 s and is large enough that 
it is accompanied by a light flash. A gas-phase breakdown is a threshold process, hence, 
it occurs when the local electric field exceeds a value characteristic of the special local 
conditions. 

A streamer develops in the wake of a positively charged tail of a primary avalanche 
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during which an accelerated electron travels through an ambient gas. If the energy of 
the accelerated electron exceeds the ionisation potential of the gas, the growth of a thin 
ionisation channel (streamer) between two electrodes results. These electrodes can be 
capacitor plates, charge distributions inside clouds, or grains. The transition between an 
avalanche process and the streamer mechanisms is not well documented so far. Raether 
(1964) suggests that the streamer mechanism dominates at an electrode distance of 5-6 cm, 
the exact value of which might change depending on local conditions. 

The onset of electric field breakdowns has been studied for charge distributions on 
capacitor plates for a variety of gases and characteristic quantities, and the breakdown 
voltage has been parametrised for easy use based on lab experiments (Eq. [I[ [2] e.g. Raizer 
1991). Different electric field strengths inside a capacitor were tested for different gases. 
For practical reasons, the parametrisation of the breakdown voltage has only been done for 
terrestrial or solar system gas mixtures or for homogeneous gases, which may not resemble 
the gas composition in an extrasolar planetary or Brown Dwarf atmosphere. 

The resulting Paschen curve, that relates the breakdown voltage with the product pd 
(Eqs. [TJ Fig. [I]), shows that the breakdown voltage of a capacitor, V t [V], depends only 
on the material of the cathodes (characterised by Townsend's second ionisation coefficient 
7), the gas between the capacitor plates (characterised by the coefficients A [cm^ 1 torr -1 ] 
and B [Vcm -1 torr -1 ]), and the product between the gas pressure and the capacitor plate 
separation, pd [cm torr], 

B-pd , „ , A 

V t = - — -V— with C = In 



C + Hpd) ~ ln(i + l) 

Et B 



(2) 



p C + \n(pd) 

Townsend's second ionisation coefficient 7 is the effective secondary emission coefficient of 
the cathode and can be calculated as the ratio between the number of electrons ejected per 
incident ion. We will refrain from doing so and apply numbers given in the literature. 7 
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varies between 10 -9 for organic cathode materials and 10 -2 for nickel and copper cathodes 
(Table 4.10 in Raizer 1991). If 7 = 1 there is one ionisation event per hit. However, the 
breakdown voltage becomes independent of the cathode material when a spark breakdown 
or a streamer starts to develop. This is evident from the fact that lightning in Earth is 
not triggered by electron emission of the cloud 'cathodes', the negatively charged part of a 
cloud. 

Qualitatively, the Paschen curve predicts that the breakdown voltage decreases as the 
electrode gap decreases because the corresponding electric field increases. The Paschen 
curve proved to be accurate for large gaps and low pressures but it is observed to fail at 
extremely low and high p<i-values. We note that experiments show that the Paschen curve 
is not necessarily accurate in describing breakdowns between electrode spaces of less than 
15/im. Go & Pohlman (2010) have therefore derived a model for a modified Paschen curve 
for breakdown in micro-scale gaps. In their model, the breakdown voltage does not have a 
singularity (+00) for small gap sizes but instead decreases further (see their Fig.l). The 
observation that pd = const does not hold for small gaps, and p and d become independent 
variables. For the small-scale cases considered, i.e. the distance between two charge centres 
in the form of grains, the gap size is just large enough to allow us to apply the classical 
Paschen curve formalism (Fig. [2]). 

At some gap size, the classical Paschen curve has a minimum, the Stoletow point, 
which is calculated from dV t /d{pd) = 0, resulting in 

pd=-H l - + l). (3) 

A 7 

pd is a local constant for a given gas composition (expressed by the constants A, B) and for 
a given capacitor surface material (expressed by the Townsend's second ionisation coefficient 
7) as long as the gap d is large enough. The product pd will therefore be the same for any 
planet or Brown Dwarf given they have the same gas composition. The capacitor plate 
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distance, d, which in our small-scale case is the distance between passing grains, changes 
with the local gas pressure, which in our case is a proxy for the vertical position in the 
atmosphere (Fig. [2|. 

Figure [I] shows V t (pd) with the Stoletow point of these gas-chemistry dependent 
Paschen curves indicated by octaeder symbols. Brown dwarfs would be best represented by 
the Paschen curve for H2 (brown dashed line) and giant gas planets will very likely exhibit 
a more complex gas-phase composition because they form from a protoplanetary disk at 
different distances from the host stars therein. The minimum breakdown voltage Vt,min [V] 
follows by inserting Eq. [3] into Eq. [T] (left), 

VWG^) = ^m(- + l). (4) 

A 7 

Thus, for each gas composition, there is a minimum voltage required to achieve electric 
breakdown in the gas. Only then can a streamer or subsequent lightning develop. It may 
however be more useful to express this condition in terms of the electric field strength E. 
With E — \E\ — I — d(j)(r)/dr\, and 0(r) the electric potential and voltage V = A<p(r) 
being the potential difference (A) between two points (hence, the units are the same), it 
follows E = V/r for a spherical charge distribution at a radial distance r from the charge 
Q (E = V/d for a capacitor). Inserting Eq. [3] into Eq. |2j the breakdown or ignition field 
strength follows with E = \ — d<p/dr\ from the following relation 

^ = B =► E tMn = Bp. (5) 
V 

Hence, (£' tjm i n /p)=const is a universal constant of a gas in the capacitor or the gas exposed 
to the electric field of some charge distribution. -Et.mm is the minimum electric field 
strength needed to achieve breakdown. Raizer (1991) and Sentman (2004) have published 
(-E't,min/p)- v alues for various gases and gas mixtures potentially representative of the solar 
system planetary atmospheres. These values are based on either experiments or calculations 
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that take into account details of the interaction processes of the electron with the ambient 
gas. For example, Sentman (2004) deals with gas mixtures by Boltzman modelling of the 
electron distribution function for a given gas composition (added to Table [TJ. Hence, details 
on the acceleration process, scattering and energy losses/gains are 'hidden' in the measured 
material constant. A detailed review of a kinetic approach can be found in Roussel-Dupre 
et al. (2008) . We summarise the values in Table [l] as some of the original literature is hard 
to come by. For more details on composition dependencies consult Chapter 6 in Meek & 
Craggs (1978). 

In summary, once the gas composition is defined for a given Brown Dwarf or gas giant 
planet, the pressure variation with height defines the variation of the minimum electric field 
needed for breakdown in the atmosphere which we will investigate in more detail below. 

3. Results 

The breakdown quantities introduced in the previous section will now be evaluated for 
a set of atmosphere structures that were previously calculated using the Drift-Phoenix 
model atmosphere code by Witte, Helling & Hauschildt (2009). These model atmospheres 
provide the required input quantities for the following calculations which are the local 
gas pressure, p ga s; and grain sizes, a, for each height in the atmosphere. Two models are 
studied in detail: T c g = 1600K, log(g)=5.0 as the example for a Brown Dwarf atmosphere, 
and T e s = 1600K, log(g)=3.0 as the example for a non-irradiated giant gas planet. Both 
atmospheres have an initially solar metallicity. Later in the paper, we compare these results 
to models of different effective temperature T cS = 1600, 1900, 2100, 2500K and 2800K 
(Fig. [8]), and the initial metallicity is varied between [M/H]= —5 ... with [M/H]=0 
indicating the solar values (Fig. [9]). The comparison amongst these models allows us 
to identify trends in the electric breakdown quantities in substellar and low-metallicity 
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Table V. Parametrisation of electric breakdown field according to Paschen curves: Material 
constants A & B taken from Raizer (1991) and from Sentman (2004), pd and V tim i n calculated 
from Eqs. |3J g 7 = 0.001 



A 

[l/(cmTorr)] 
(Raizer 1991) 



B 

[V/(cmTorr)] 
(Raizer 1991) 



Papain 

[cm Torr] 
(this paper) 



[V] 

(this paper) 



He 
H 2 
N 2 
Air 
C0 2 
H 2 



3 
5 

12 
15 
20 
13 



34 
130 
342 
365 
466 
290 



4.181 
2.509 
1.045 
0.8363 
0.6272 
0.9649 



Sentman (2004) 



A 

[l/(cmTorr)] 



E/p=B 
[V/ (cm Torr)] 



pdmin 

[cm Torr] 



Venus (C0 2 /N 2 =96.5/3.5) 
Earth (N 2 /0 2 /Ar=78/21/l) 
Mars (C0 2 /N 2 /Ar=95.5/2.8/1.7) 
Jupiter (H 2 /He/CH 4 =89/10.9/0.1) 
Saturn (H 2 /He/CH 4 =96.3/3.6/0.1) 
Titan & Tritan (N 2 /CH 4 =95/5) 
Uranus (H 2 /He/CH 4 =82. 5/15.2/2.3) 
Neptune (H 2 /He/CH 4 =80/18.5/1.5) 



7.27 
7.44 
7.23 
6.19 
7.46 
8.80 
6.47 
622 



180 
243 
178 
143 
156 
274 
138 
135 



2.58 
2.53 
2.60 
3.06 
2.52 
2.13 
2.90 
3.02 
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atmospheres. 

3.1. Classical breakdown quantities 

Breakdown distance d: Assuming that a charge separation has built up in a Brown 
Dwarf or planetary atmosphere, how far can two imaginary capacitor plates be apart that 
a breakdown of the atmospheric gas will occur for an established electric field? 

This is the distance of two capacitor plates, or two otherwise charge-carrying surfaces, 
that can build up an electric field between them, as described by the Paschen minimum 
where p gas • d = 2.509 torr cm (=5007 dyn/cm 2 ) for an H 2 gas (Eq. |3j Fig. [I]). This is 
the distance where the electrical breakdown would be easiest according to Eq. |4j Hence, 
for p gas ■ d below the Paschen minimum, the breakdown voltage increases sharply in the 
classical case. However, Go & Pohlman (2010) show that lower breakdown voltages should 
be expected at smaller distances, justifying our use of the classical Paschen curve as a first 
step to study breakdown characteristics. p gas • d is a constant for a given A [cm^ 1 torr -1 ] 
and 7, hence, the same for all atmospheres of the same chemical composition with respect 
to the most abundant gas-phase species. This allows us to calculate the breakdown distance 
d for a known p gas -structure of a model atmosphere. 

Figure [2] shows the breakdown distance as function of the local gas pressure calculated 
for the pressure structure of a Brown Dwarf (T e g- = 1600K, log(g)=5.0, solar element 
abundances, brown line) and a giant gas planet (T c g- = 1600K, log(g)=3.0, solar element 
abundances, red lines) assuming that molecular hydrogen (A[H 2 ]=5 cm _1 Torr _1 , Raizer 
1991) is the dominating gas species to be ionised by electrons accelerated in the electric field. 
This assumption is appropriate as H 2 is the most abundant molecule in the atmosphere 
models applied here. The assumption breaks down if H 2 ionises in the hot and dense inner 
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layers of the atmosphere where clouds can not be thermally stable, and where atomic 
hydrogen is the dominating gas component. The Townsend's 2nd ionisation coefficient 
7 = 0.001 is used, representing a conducting cathode material which is strictly correct only 
for T >1800K from which iron dominates the grain material composition in the cloud. 
As mentioned earlier, it is only relevant to consider the material dependence of 7 for the 
avalanche process, not so for the streamer. 

Figure [2] shows that the breakdown distance varies over the atmosphere pressure 
ranges from ~ 5/zm in the densest regions to 10 5 km in the upper, low pressure regions 
(dashed lines). This suggests that discharges could appear spatially extended and diffuse 
in low-pressure regions, while numerous discharge events remain small, locally confined 
phenomena in high-pressure regions. It also suggests that different electric discharge 
processes dominate at different heights: local corona (point discharges) and small-scale 
sparks at the bottom of the cloud, and glow discharge and large-scale sparks at the top. 

In a quasi-static environment, these different processes could lead to a neutralisation 
of large-scale charge built-up and consequently diminish the possibility for a large-scale 
lightning. However, hydrodynamic atmosphere motions are driven by convection and by 
gravity waves (Freytag et al. 2010) on Brown Dwarfs which are also rapid rotators (Scholz 
et al. 2011). This suggests that large-scale charge separation is likely to be sustained by 
convection-like processes comparable to Earth, for example like in the Hadley cells at the 
Earth's equator. 

Following the observation in Raizer (1991) that for breakdown distances > 5 . . . 6cm 
the streamer mechanism will dominate and sparks can develop, we distinguish a large-scale 
(d> 6cm) and a small-scale regime (d< 6cm) in Fig [2] 

It is interesting to compare the maximum spacial separation that two charge 
distributions whilst still permitting a discharge, to the average size of the cloud particles 
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and their mean particle distance. For simplicity, we consider the altitude dependent 
mean grain size (a) instead of separate grain size distributions for each altitude as shown 
in Helling, Woitke & Thi (2008). The mean grain size, (a), is always smaller than the 
maximum distance, d, over the whole geometrical extension of the cloud. The mean 

— 1/3 

distance between the cloud particles scales with the cloud particle number density as n d . 
Consequently, it will be largest where the number density is smallest. The 'knee' at the 
high-pressure end of each long-dashed curve in Fig. [2] corresponds to the maximum cloud 



particle number density (compare the upper panel in Fig. 10). Generally, the cloud particles 
have a mean separation that is ~ 3.5 orders of magnitude smaller than the maximum 
break-through distance in a giant gas planet cloud (~ 2.5 for the Brown Dwarf cloud). This 
comparison also shows that the cloud particles have considerable separation in the upper 
cloud region suggesting that collective phenomena such as large-scale charge distributions 
are to be expected rather than inter-particle discharges such as in the high-pressure part of 
the cloud. However, if the breakdown distance is smaller than the maximum value d (e.g. 
in a turbulent flow), small-scale discharge processes may also occur in the low-pressure part 
of an atmosphere. 

Note, the mean particle size and the cloud particle number density are determined by 
the processes of dust and cloud formation (nucleation, growth/evaporation, drift, convective 
mixing, element depletion) taken into account in the Drift-Phoenix model atmospheres. 
Coagulation, for instance, could change the grain size distribution, but it is not clear a 
priori how it changes since a variety of constructive and destructive collisional processes 
occur (Guttler et al. 2010). Furthermore, turbulence (Helling et al. 2001, 2004) and gravity 
waves (Freytag et al. 2010) cause local, hydrodynamic density increases, neither of which is 
taken into account in the present paper. 

Depending on the ionising gas, these global values of d have uncertainties of one 
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order of magnitude as we demonstrate for different gas-phase compositions in Fig. [3] The 
breakdown distances inside the cloud regions (solid lines, Fig. [2J are somewhat different 
in a Brown Dwarf mineral cloud and in a giant planet's mineral cloud because the Brown 
Dwarf clouds exist at higher pressures. The exact values of discharge quantities, such as 
the maximum separation of charge distributions before breakdown but also the required 
number of charges, will therefore differ between objects of different effective temperature, 
surface gravity and metallicity (see Figs. [8j [9]). 

Minimum electric field for breakdown E t<min : The minimum electric field strength, 
needed for a gas breakdown scales with the local gas pressure (Eq. [5]), are shown in Fig. |1J 
This critical breakdown field strength corresponds to the minimum of the Paschen curve for 
individual ionised gases (H 2 , He, air, H 2 0, CO2; Table 4.1 in Raizer 1991) represented by 
individual parameters B [V cm -1 torr -1 ], and a given Townsend's 2nd ionisation coefficient 
7 = 0.001. The results are shown for a sample model atmosphere (T e ff=1600K, log(g)=5, 
initial solar composition) for which we vary the parameter B to show the dependence on gas 
composition. Figure [4] shows results for pure gases (data from Raizer 1991), and Fig. [5] for 
mixed gases of the solar system planets (data from Sentman 2004). Over-plotted in Fig. [4] 
are the values given in Yair et al. (1995: magenta - Jupiter cloud with drops and ice, cyan 
- Earth values; full symbol - Paschen formula, open symbols - measurements). This shows 
that the values are comparable. The vertical lines in Fig. [4] indicate the large-scale regime 
(low pressure) and the small-scale regime (high pressure) for pressures corresponding to 
Fig-U 

The minimum critical breakdown field strength varies between 10 7 V/cm at high 
atmospheric pressures and 10 -7 V/cm in low pressure regions where p gas = 10 -12 bar. The 
breakdown fields on Earth (at sea level) and in Jupiter are ~ 10 4 . . . 10 6 V/cm which 
correspond to the values at 1 bar in Figs. [4] and [5} The dependence on the local gas 



-19- 



composition introduces an uncertainty of one order of magnitude as demonstrated in Fig. |4j 

Figure [5] shows a close-up of Fig. [4] showing the minimum breakdown field using 
the data for the different solar system planets in comparison with the two most different 
examples of gas phase composition (CO2 and He gas). The results for the solar system 
planets sit between a pure H 2 (upper limit) and a pure H 2 atmosphere (lower limit). 

We note further that values measured for the onset of coronal discharges on a treeless 
plot of grass are about 3 orders of magnitude lower (40 V/cm; MacGorman Sz Rust 1998) 
than the respective value at 1 bar in Figs. |4j [5j MacGorman & Rust (1998) also point out 
that point discharges from ice crystals which are dominated by their surface (rather than 
bulk) conductivity, seem likely in thunderstorms for a field strength > 400 kV/m (4 • 10 3 
V/cm) which is still ~ 2 orders of magnitude lower than the classical values in Figs. [4] and 

El 

Therefore, our results suggest that electrical breakdown can readily occur in cloudy 
atmospheres of Brown Dwarfs and giant gas planets. 

3.2. How many charges are needed for an electric field breakdown? 

Laboratory measurements of grain charges suggest that single grains can carry 
10 3 — 10 5 e/grain. This is supported by experiments on photoelectrical charging of dust in 
a vacuum (Sickafoose et al. 2000; metals and glass), grain levitation (Fortov et al. 2001), 
measurements of grauple on Earth (Lamb & Verlinde 2011), and wind-driven entrainment 
of Martian dust (Merrison et al. 2012; silicates and iron oxide). Somewhat higher charge 
loads are suggested from dust-dust collision experiments (e.g. contact electrification and 
tribo-electrification) with silicate granules by Poppe, Blum & Henning (2000) resulting in 
10 -5 C/m -2 (i.e., ~ 5 • 10 6 e/grain for a=1.2/mi), and by volcanic ash experiments involving 
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dust charging by fractoemisson by James et al. (2000) resulting in 10 -8 . . . 10 -6 ' 5 C/m -2 
(~ 10 5 . . . 10 7 ' 5 e/cm 2 ). How do these experimental numbers compare to the number of 
charges needed to initiate a field breakdown in a substellar atmosphere? 

The minimum electric field strength that needs to be overcome in order to produce a 

upper 



field breakdown in an astrophysical gas with a low degree of ionisation (see Fig. 10 
panel, solid lines) varies only by one order of magnitude between gases of different molecular 
composition for a given gas pressure. Furthermore, the breakdown voltage measured above 
thunderclouds on Earth is up to two order of magnitude lower than the classical breakdown 



values that were discussed here so far (see also the discussion at the end of Sect. 3.1). 
Keeping these uncertainties in mind, we continue with studying how many charges are 
needed to initiate a field breakdown in a gas of a given composition for the minimum 
breakdown voltage K, m in ■ 

The assumption of a spherically symmetric charge distribution (Eq. [6]) allows a first 
order-of-magnitude study, although neither clouds nor cloud particles are truly spherical. 
Equation [6] is a general expression which allows us to distinguish between these two scale 
regimes of clouds and grains later on. With the radial distance from a charge distribution 
of Qcrit charges set to the minimum distance for breakdown, d, it follows that 

Q crit = 47re dV tMn [C] Q crit oc d oc — ^r- (6) 

1 electron carries only the tiny fraction of a Coulomb: 1C = 6.24150965 10 18 e. The surface 
charge density for a charge distribution of radius r (grain radius or cloud radius) is given 
by (7 = Q/(4irr 2 ). Thus, the critical surface charge density for breakdown is 

<xcnt = e -^^ = ^E„(^) 2 [C/cm 2 ] (7) 

(eo = 8.85 10 _12 Fm _1 (r]) - vacuum dielectric constant or electric permittivity of free space). 



2 e is in SI unites Fm" 1 with F=As/V (A - Ampere) = C/V (C-Coulomb) = C 2 /(Nm) 
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Note that the proportionality at the right hand side of Eq|6] only holds if the radial distance 
from the charge distribution is set to the minimum distance for breakdown according to the 
minimum of the Paschen curve (Eq. [3]). Note further that the pressure needs to be in SI 
units [Pa] for Q C nt to be in units of elementary charges on the right hand side of Eq. [6] 

The charge carrying surface, 4irr 2 , of our charge distribution can be the effective 
surface of a cloud or the surface of a single grain (r = a). The total charge carrying surface 
of the cloud, the effective cloud surface, will be provided by the total surface of the particles 
that form the cloud, rather than being related to the geometrical cloud radius. On large 
vertical scales, the cloud particles will settle depending on their masses and sizes. Our 
stationary cloud formation model suggests a large mean grain size at the lower edge of 
the clouds and a small mean grain size in the cloud's top region (Woitke & Helling 2004; 
Witte, Helling & Hauschildt 2009). Smaller particles will remain suspended for longer in 
the atmosphere than the more positively charged bigger grains. On small scales, the charge 
separation needed to build up an electric field is influenced by non-spherical grain shapes 
and chemical surface inhomogeneities, neither of which is included in our considerations 
yet. For spheroidal grain growth, we refer to Diver & Clark (1997) and Stark et al. (2006). 

Our cloud formation model (Woitke & Helling 2003, 2004; Helling & Woitke 2006; 
Helling, Woitke & Thi 2008) allows us to link numerically derived cloud properties, like 
grain sizes, grain number densities, and dust surface, to critical breakdown quantities as 
evaluated in the previous sections. The solution of the our cloud model in the frame of the 
radiative transfer code Phoenix (Hauschildt & Baron 1999) allows to link dust properties 
and the critical breakdown quantities to the atmospheric temperature and pressure scale 
of a Brown Dwarf or giant gas planet as result of the Drift-Phoenix model atmosphere 
simulations (Dehn 2007, Helling et al. 2007a,b; Witte, Helling & Hauschildt 2009). 

The moment equations of our kinetic dust formation model (Gail & Sedlmayr 1988, 
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Dominik et al. 1993, Woitke & Helling 2003) allow us to calculate the total dust surface 
available per unit volume (^4ust)' anc ^ a ^ so the mean surface of a grain (<v4 dust >) depending 
on atmospheric height (i.e. depending on the local temperature and pressure), 

AX st = v / 36VL 2 [cm 2 /cm 3 ] (8) 
<A dust > = ^3&KL 2 /L Q = A^Jn d [cm 2 ]. (9) 

Lq is the zeroth and L 2 the second moment of the grain size distribution function, 
respectively (e.g. Helling et al. 2001). Equations [8] and [9] allow us to relate the total number 
of charges needed for breakdown to the total grain surface per unit volume (upper panel, 
Fig. [7]) and to the numbers of grains (lower panel, Fig. [7]) available as follows: 

a x A d ^. t [e~/cm 3 ] — number of charges per total dust surface per gas volume (10) 
a x dust [e - /grain] — number of charges per grain. (11) 

Note that the total dust surface, i.e. the total cloud particle surface, is related to the charge 
carrying surface inside a cloud. 

The number of charges required to establish an electric field large enough for a 
breakdown to occur and subsequently to ionise the ambient gas, depends on the local gas 
pressure. More charges are needed in the low density (pressure) regime of an atmosphere 
as for example Fig. [6] demonstrates because of the large maximum breakdown distance 
(compare Figs. |2j [3]). The lowest density layers would require a total of about 10 20 charges 
to allow an electron avalanche with a subsequent streamer to develop. The actual number 
depends to a certain extent on the dominating gas species that is being ionised as in the 
example for H 2 and CO2 in Fig. |6} The lower panel of Fig. [6] presents this results in (SI) 
units of critical charge density, <7 cr i t [C/cm 2 ], assuming that the charges would be located 
on the grains of the size of the local mean grain size, hence r =<a(T, p gas )>. Artificially 
increasing the radius of the charge carrying surface to r = 100x<a(T, p gas )> (long dashed 
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red line) verifies the anticipated result of a decreasing critical charge density with increasing 
surface area that would allow an electric field breakdown. 

Figure [7] (upper panel) shows the critical number of charges, a cvl t x A^ st [e~/cm 3 ], as 
distributed over the total dust surface per cm 3 of atmospheric gas depending on altitude 
(represented by the local atmospheric gas pressure). The dust surface is provided by a 
large number of large dust grains at higher pressure and by small number of smaller grains 
at lower pressures at the cloud top. If the same total dust surface would be provided by 
grains of a constant size of 10 _3 cm (brown long-dash short dash line), the critical number 
of charges per cm 3 would decrease by more than 5 orders of magnitude. 

shows the critical number of charges per grain, cx cr it x -^j L 
[e~ /grain], which takes into account the total dust surface and the number of grains present 
throughout the cloud. Each grain has to carry considerably less charge at the bottom of 
the cloud compared to the cloud top as here more grains are available. Due to gravitational 
settling, most of the grains are contained in the deeper cloud where the number of grains is 
much higher than at the cloud top. 

Our results suggest that many charges per grain are required to initiate a breakdown 
discharge in the upper part of the clouds in Brown Dwarfs and giant gas planets. Such a 
height-dependent charge distribution results because we apply our results for the maximum 
breakdown distance, d, which is rather large at the cloud top as demonstrated in Figs. [2j [3] 
The breakdown distance can, however, decrease by e.g. turbulent or other non-local wind 
processes which are not taken into account here. A decreasing number of charges needed 
to initiate an electric field breakdown would then result. Additionally, the charge density, 
a, depends on the radius of the charge carrying surface, r. Therefore, we again assess our 
results for r larger than <a>, but now by using a constant grain size of 10~ 3 cm (brown 
long-dash short dash line) throughout the cloud. The comparison with the experiments 
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with mono-disperse grain ensembles as performed by Sickafoose et al. (2000; green star in 
Fig. [7]), Fortov et al. (2001; magenta vertical bar in Fig. [7]), Merrison et al. (2012; blue box 
in Fig. [7]), Poppe, Blume & Henning (2000; blue asterisk in Fig. [7]), James et al (2000; blue 
asterisk in Fig. [7]), and values for grauple on Earth from Lamb & Verlinde (2011; orange 
vertical bar in Fig. [7]) suggest that our results for mineral clouds are not unreasonable, 
especially if we allow that cloud particle sizes can be larger than the mean grain size from 
our dust model. Such a deviation from the mean grain size is to be expected as our study 
of the grain size distribution in Helling, Woitke & Thi (2008) demonstrates. Furthermore, 
critical breakdown field strengths applied here are upper values as experiment on Earth 
thunderstorms suggest and can, hence, be expected to be 1-2 orders of magnitude lower 
with the corresponding decreasing effect on the necessary number of charges. The chemical 
richness of the mineral clouds in the atmosphere of Brown Dwarfs and giant gas planets 
may introduce further uncertainties compared to the laboratory experiments cited above. 
The actual material composition may influence the charge distribution on the grain surface 
leading, for example, to further polarisation effects of the cloud particles. The question of 
how many charges a dust grain can carry before it is destroyed by the electrostatic stress 
will be addressed in a forthcoming paper. 

Trends with T c g and metallicity [M/H]: The quantities that determine if the streamer 
mechanism occurs and with which efficiency it does so, are determined by the local gas 
pressure because the gas pressure determines the mean free path of the seed electrons inside 
the gas of Brown Dwarfs and planets in the electric field. Global quantities that determine 
an atmosphere's pressure structure are the effective temperature, T eff , the metallicity, 
[M/H], and the surface gravity, log(g), of the object, in the case of non-irradiated objects. 
The influence of the surface gravity has already been shown in Figs. [6] and [7j We use the 
surface gravity to distinguish between Brown Dwarfs (log(g)=5.0) and giant gas planets 
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(log(g)=3.0). 

Figures [I] and [9] show the influence of T e fr and [M/H], respectively, for a small 
number of Drift-Phoenix models. The metallicity [M/H] is given relative to the solar 
abundances, hence, [M/H]=— 3.0 stands for the initial element abundances being decreased 
by 10~ 3 compared to the solar values. The lowest panels in both figures demonstrate 
how the Drift-Phoenix temperature-pressure profiles change with changing metallicity. 
The local pressure mainly increases for a given local temperature for an increasing T e jf 
(Figures [8j), while the upper atmospheric parts of the (T gas , p gas )-structures stay relatively 
comparable with decreasing metallicity for given T e ff and log(g). The inner atmospheric 
regions, however, show an increasing gas pressure with decreasing metallicity for a given 
gas temperature. 

The critical charge number density needed to initiate the occurrence of streamers 
that subsequently might lead to lightning, does change by up to 4 orders of magnitude 
for a given gas pressure amongst the models of varying effective temperature that form 
dust clouds. Note that the T e g- = 2800K-model in Fig. [8] is too hot for dust condensation 
processes to be efficient. The metallicity of the atmospheric gas has a larger impact on 
the local temperature. Figure [9] shows that the critical charge number density can differ 
in the low-metallicity cases much more than for varying T e ff. Note that the dust surface 
available to host the charges does vary amongst the models, with the largest differences 
occurring again for changing metallicities in the lowest metallicity cases (dashed/dotted 
line) compared to the solar case (solid line). All models, however, exhibit a comparable 
total dust surface in the inner, denser regions at the bottom of the cloud. Witte et al. 
(2009) demonstrated that low-metallicity atmospheres can produce larger dust grains at 
the inner edge of the cloud than in the solar case if [M/H]>=-4.0. 
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3.3. Number of electrons vs. number of grains 

So far, we have studied the breakdown conditions utilising classical arguments and 
laboratory measurements assuming that the gaseous atmosphere would have enough free 
electrons available for starting the electric breakdown process. But, is the number of 
thermal electrons large enough in Brown Dwarf and giant gas planet atmospheres to 
provide the seed electrons that are needed to start an electron avalanche in the electric 
field, for example between two passing grains? The investigation of this question allows us 
to estimate the efficiency with which streamer discharges could occur inside Brown Dwarf 
and giant gas planet atmospheres. Streamer discharges are the starting point for large-scale 
phenomena like lightning that may be followed by a sprite into the less dense atmosphere 
above. 

We compare the degree of thermal ionisation, Xe,themai = ^e,thermai/^H (solid lines), with 



the number of dust grains per hydrogen, n^/nn (dashed lines) in Fig. 10 (top panel), and we 



use these two quantities to estimate an efficiency of streamer initiation in a mineral cloud 



by calculating the ratio n ethermal /n d (Fig. 10, bottom panel). Hence, this efficiency would 



be the number of streamers per grain pair if each of the thermal electrons in the electric 
field of two charged grains would initiate an electron avalanche. In the previous section, we 
have evaluated the conditions under which such a field breakdown can occur. 



Figure 10 shows that the efficiency, n e ,thermai/^d) varies with height as at some height 
the dust number density can have a local maximum while the number of thermal electrons 
is a rather smooth function of gas pressure which is linked to the gas temperature by the 
underlying atmosphere model. Hence, the number of electrons available for each colliding 



pair of charge-carrying grains decreases as shown in the lower panel of Fig. 10 Here, more 
than one electron would be available for initiating the avalanche process for all pressures 
above the solid black, horizontal lines (^ e ,thermai — n d), an d less than one electron would 
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be available below this line. In both cases, the Brown Dwarf atmosphere (brown) and 
the giant gas planet's atmosphere (red), the efficiency of streamer initiation would vary 
between 1 and ~ 1CT 5 . n ei thermai/^d = 1 means that one streamer per grain pair sets off, 
and n Ci thermai < n<i means that not every grain pair can initiate a streamer event. To 
generalise, we can distinguish two regimes regarding the streamer efficiency in a cloud 
forming atmosphere, the occurrence of which can be less homogeneous if the whole grain 
size distributions for each atmospheric layer are taken into account, 

n e ,thcrmai > n d electron dominated n e ,thcrmai/ n d > 1 

^thermal < ™ d dust dominated ™e,thermal/™d < 1 • • • 10" 5 . 

The consequence is that the discharge process in cloudy atmospheres will be determined 
by two populations, the population of the cloud particles and the population of free electrons, 
respectively. Both populations need to be abundant enough to allow discharge processes to 
develop into streamers as the pre-condition for the occurrence of lightning. 



3.4. An estimate of electron enrichment of the gas phase by dust-dust induced 

streamer events in mineral clouds 

The empirical formula 

a = Apexp(-^) a(E min ) = ^ (12) 

hi e 

for Townsend's ionisation coefficient, a [cm -1 ], enables us to provide first estimates of the 
number of electrons created when a electron avalanche starts that might develop into a 
streamer. Streamers are discussed to be a likely possibility to initiate lightning (see Sect. 
5.4. in MacGourman & Rust 1998). a is the number of ionisation events performed by 
an electron in a 1cm path along an electric field. This interpretation is only true if we 
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assume that the electron undergoes ionisation collisions only for high values of E/p and 
moderate electron energies (Raizer 1991). The constants A and B are the same constant 
given in Table [T] for approximating experimental curves, ad would then be an estimate for 
the number of electrons produced if the electron avalanche would stretch along the whole 
critical breakdown distance d (compare Fig. |2|. The total number of non-thermal electrons 
per unit time, Aff^fs -1 ], produced by all possible seed electrons that lead to an electron 
avalanche, can be approximated by 

iV e aval « fes xadx^ n (z) (13) 

^e,thermal i dd / \ /i »\ 

sa -j- x ad x uZi(z), (14) 

i/z ■ tld 

with t'coiil' 2 ) the frequency of dust-dust collisions at each height z in the atmosphere as 



derived in Eq.17 in Helling, Jardine & Mokler (2011). Following from Sect. 3.3 we define 
fcs = thermal/ (1/2 ' n d) which provides an estimate of the efficiency with which an 
avalanche could occur. f e s > 1 suggests a superposition of avalanches that are each started 
by one gas-phase electron. 

An electron avalanche can turn into a self-propagating ionisation front by the 
subsequent development of a streamer in its positively charged wake. This instability can 
trigger much larger plasma-channel structures like sprites above thunderclouds or lightning 
bolts below and inside a cloud which we are not considering for our estimate here. The 
elementary streamer alone, however, results in a number of local non-thermal electrons much 



larger than predicted by the classical Townsend breakdown by electron avalanche (Eq. 12). 
Numerical experiments suggest the production of 10 12 . . . 10 13 electrons per streamer event 
(e.g. Dowds et al. 1993). Multiple streamers develop if more than one electron is 'seen' 
by a pair of charged grains or two otherwise oppositely charged charge distributions. This 
efficiency depends on the number of free electrons that are already present in the gas before 



a streamer-triggering electron avalanche develops and has been derived in Sect. 3.3 For 
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a streamer breakdown, the number of non-thermal electrons produced during a collisional 
approach of two dust particles that occur with the frequency v^y^z), 

iVf— = ^xl0 12 x^). (15) 
1/2 • n d 

Given the different nature of an electron avalanche and a streamer, the electron 
production occurs with a different efficiency. The result is that jV e streamer > iV aval (Fig. m}. 



Our results in Fig. 11 suggest further that the largest impact of the streamer electrons on 
the local chemistry could be expected in the inner, denser cloud regions in Brown Dwarfs 
and giant gas planets. 



4. Conclusions 

Discharge process require a) charge separation on cloud particles and b) a charge 
separation over large distances to allow large-scale discharges such as lightning or sprites 
to occur. The laboratory experiments (Poppe, Blum & Henning 2000, James et al. 2000, 
Fortov et al. 2001, Lacks & Levandovsky 2007, Forward et al. 2009, Merrison et al. 
2012) for a variety of minerals (including Martian and volcanic minerals), for graupel 
on Earth (Lamb & Verlinde 2011), as well as our own investigations (Helling, Jardine, 
Mokler 2011) suggest that mineral cloud particles will be charged. Mineral cloud particles 
settle gravitationally in Brown Dwarfs and giant gas planets such that a large scale charge 
separation can occur inside these clouds. Convection serves as an additional large-scale 
charge separation mechanism. If small grains are indeed more negatively charged than 
large grains as suggested by experiments (Hatakeyama & Uchikawa 1952, Kikuchi & 
Endoh 1982, James et al. 2000, Merrison et al. 2012), spatial charge separation will occur 
on a large scale by gravitational settling but also on small scales by particle size effects 
alone. Our investigations of the minimum distance that two charge distributions need to 
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have for a electric field breakdown to occur (Fig. [2]) suggest that small-scale discharges 
should be expected in high-pressure regions of an atmospheric cloud. Large-scale discharge 
processes should occur in the upper, low-pressure part of the cloud. This is consistent with 
observations from Earth where, for example, sprites develop into the low-pressure part of 
the atmosphere above a thunder-cloud event. 

The critical electric field strength to be overcome for a field breakdown varies between 
10 7 V/cm at high atmospheric pressures and 10 -7 V/cm in the upper atmosphere where the 
gas pressure is low. The critical number of charges per total dust surface per cm 3 change 
from 10 23 e/cm 3 in the inner atmosphere to < 10 5 e/cm 3 in the low pressure atmosphere. 
These numbers compare well with experimental values at the lbar-pressure level. The 
number of charges needed to initiate the occurrence of streamers that subsequently might 
lead to lightning, does not change substantially amongst the atmosphere models of varying 
effective temperature that form dust clouds. The critical number of charges has the same 
trend in the low-metallicity cases as for varying effective temperatures. 

Generally, a population of charged cloud particles and enough free electrons are needed 
to allow transient luminous events to occur in extraterrestrial, cloudy atmospheres. The 
ratio of both determines the efficiency with which discharge events may occur. The rate of 
electron enrichment is highest in the high-pressure part of the cloud as here, according to 
our model, a larger charge- carrying surface is available. This is where we expect small-scale 
discharges to dominate. 
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Fig. 1. — The Paschen curve shows the dependence of the breakdown voltage, Vt [V] on the product 
of gas pressure, p, and the minimum breakdown distance, d, for four molecular gases composed of 
H2, H2O, CO2, and of air. The classical Paschen curve depicted here, has a minimum (indicated 
by the symbol; Stoletow point) which we use to determine the minimum electric breakdown field 
in Brown Dwarf and exoplanetary atmospheres in Figs. [2l [3] 
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Fig. 2. — The minimum separation of charges of two capacitor plates (representing our charge 
carrying surfaces like on clouds or on grains) required for discharge (solid and short-dashed in the 
upper part of the plot). The distance d corresponds to the Paschen minimum where p gas ■ d = 2.509 
torr cm (=5007 dyn/cm 2 cm) for an H2 gas. Note, an artificial offset between the log(g)=5.0 
(brown, "logg=5.0") and the log(g)=3.0 (red, "logg=3.0") was applied; all curves should lie on top 
of each other without this offset. p g3bS ■ d is a constant for a given A [cm -1 torr -1 ] and 7, hence, 
the same for all atmospheres of the same chemical composition. In this plot, an H2-gas is assumed 
for which A[H2]=5 cm -1 Torr -1 (Raizer 1991) and Townsend's 2nd ionisation coefficient 7 = 0.001 
is used. The solid lines show the pressure interval where clouds form and the dashed lines show 
the whole atmospheric pressure range modelled. Raizer (1991) notes that discharge sparks set in 
if the plate distance is >6cm (horizontal dotted line) which helps us here to distinguish between 
large-scale and small-scale discharges. We overplot two more length scales for comparison: The 
dust-dust mean free path n d (long dashed lines) and the mean grain sizes, (a), at the bottom of 
the plot. 



-37- 



10 - 



O 



-i 1 1 1 1 1 1 1 1 1 1 1 1 1 

^ Teff=1600K, log(g)=5.0, solar 

7= 001 
H2 

\^ - He 

' • V N - Air 

'^V H20 

•:-^ s C02 

X\ 

■v, 

\\ 

X\ 

_l I I I I I I I I I I I I 

-10 -5 

log P gM [bar] 



Fig. 3. — Dependance of breakdown distance d at the Stoletow point on the composition of the 
ionised gas. The variation in these distances for 6 cases (H2, He, air, H2O, CO2; Table 4.1 in Raizer 
1991) shows that the uncertainty resulting from atmospheric composition is a factor of 10. This 
estimate also covers differences in A depending on the pressure regime in which the measurements 
were made (high [~ 1 bar] vs. low pressure; Table 7.1 in Raizer 1991). Sentman (2004) has 
published data for mixed gases but the lack of a consistent notation makes application difficult. 
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Fig. 4. — The minimum critical breakdown electric field corresponding to the minimum of the 
Paschen curve for individual ionised gases (H2, He, air, H2O, CO2; Table 4.1 in Raizer 1991). 
Pressures higher than indicated by the vertical dotted line correspond to d < 6cm, and pressures 
lower than indicated by the vertical dotted line correspond to d > 6cm according to Fig. [3j (This 
set of dotted double lines results from the intersection interval of the dotted line in Fig. |3j) The 
symbols indicate the breakdown field for Earth (cyan) and Jupiter (magenta) as given in Yair et 
al. 1997 (Table 1; open symbols: measured, solid: with drops and ice, stars: theroretical estimate). 
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Fig. 5. — The minimum critical breakdown electric field for a Brown Dwarf atmosphere with T e g = 
1600K and log(g)=5.0 assuming that the field break down happens in different gas compositions. 
This is a zoom into Fig.|4]for CO2 (red dotted), H2O (blue long-dashed), H2 (brown dashed) and He 
(orange long dashed) with the curves assuming the solar system planetary atmosphere compositions 
over-plotted in black. The curve for Uranus and Neptune lie almost on top of each other. 
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Fig. 6. — Top: The total number of charges needed to overcome the breakdown field (Eq. [6]) 
for a H2 dominated gas in a mineral cloud in a brown dwarf atmosphere (brown) and a giant gas 
planet (red). These charges are assumed to exist on two charge-carrying surfaces separated by 
the minimum distance for discharge. The case of a C02-dominated atmosphere is over-plotted for 
the giant-gas planet (dotted line). Bottom: The critical surface charge density, cr cr = Q C r/(47rr 2 ) 
critical for an electric field breakdown. cr cr was calculated for r =<a(T, p gas )> assuming that all 
grains are of the size of the local mean grain size, except for the long-dashed line which shows cr cr 
for 100x<a>. Note that 1 C/m 2 = 6.24 • 10 14 e/cm 2 . 
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Fig. 7. — Critical number of charges needed for electric breakdown in mineral clouds inside atmo- 
spheres for objects with T e g=1600K and log(g)=3.0 (red) and log(g)=5.0 (brown). Upper panel: 
The charges on the total dust surface A^ st per cm 3 of atmospheric gas as a cr x A^ st [e/cm 3 ]. 
Lower panel: The number of charges per dust grain cr cr x (^dust/ n d) [e/grain]. The symbols 
represent the following data sources: grey rectangle - Merrison et al. (2012), magenta vertical bar 
- Fortov et al. (2001), green asterisk - Sickafoose et al. (2000), blue cross - Poppe et al. (2000) 
and James et al. (2000), orange vertical bar - values for Earth hale from Lamb & Verlinde (2011). 
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Fig. 8. — Critical charge number density, a cr [e/cm 3 ] per total dust surface needed to overcome the 
local breakdown field depending on the effective temperature of the object (top panel). The results 
are plotted for 5 solar-metalicity Brown Dwarf model atmospheres (log(g)=5.0) with T e g- = 1600K 
(solid), 1900K (dotted), 2100K (short dashed), 2500K (long dashed), 2800K (dash-dot). The lower 
panel (local gas temperature) demonstrates why the critical surface charge density a cr (p) changes 
for different effective temperatures. 
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Fig. 9. — Critical charge number density, a cr [e/cm 3 ] per total dust surface, needed to overcome 
the local breakdown field depending on the metalliciy of the object (top panel). The results are 
plotted for 5 Brown Dwarf model atmospheres (log(g)=5.0) with different metallicities [M/H]=0 
(solar, solid), -1 (dotted), -1 (short dashed), -3 (long dashed), -5 (dash-dot). The lower panel (local 
gas temperature) demonstrate why the critical surface charge density cr C r(p) changes for different 
[M/H] . 
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Fig. 10. — The efficiency of discharges in planetary atmospheres with mineral clouds. Top: 
De gree of thermal gas ionisation, Xe,themai — /nn (solid lines), and the dust number density per 
hydrogen, n^/nn (dashed lines), for a Brown Dwarf and a giant gas planet atmosphere for the 
same Drift-Phoenix models as in Fig. [6| Bottom: The ratio of the number of thermal electrons, 
^e.thermai, to the number of cloud particles, n^, provides a measure of the efficiency with which a 
streamer may be initiated in a mineral cloud. Two different regimes appear: above the black line - 
electron dominated, hence more than one electron is available per grain pair, below the black line 
- dust dominated, hence less than one electron is available per grain pair. The dotted vertical lines 
distinguished the cloud regimes of potential large-scale cloud discharges and small-scale, inter- grain 
discharges similar to Figs. [2]-[4] 
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Fig. 11. — Gas-phase enrichment by dust-dust collision triggered discharges. Electron avalanche 
(dashed lines; d = 5cm used) yields considerably less free charges than streamer events (solid lines). 
The same Brown Dwarf (brown) and a giant gas planet (red) atmosphere models are used as in 
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